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Borrowing and extending the method of images we introduce a theoretical Iramework that greatly 
simplifies analytical and numerical investigations of the escape rate in open dynamical systems. As 
an example, we explicitly derive the exact size- and position-dependent escape rate in a Markov case 
for holes of finite size. Moreover, a general relation between the transfer operators of closed and 
corresponding open systems, together with the generating function of the probability of return to the 
hole is derived. This relation is then used to compute the small hole asymptotic behavior, in terms 
of readily calculable quantities. As an example we derive logarithmic corrections in the second order 
term. Being valid for Markov systems, our framework can find application in information theory, 
network theory, quantum Weyl law and via Ulam's method can be used as an approximation method 
in more general dynamical systems. 

PACS numbers: 05.45.Ac, 02.50.Ga. 



The study of the statistical properties emerging from chaotic dynamics typically deals with equilibrium quantities 
and related convergence issues, asymptotic in time. On the other hand, there is growing interest in applications 
of dynamical systems theory to problems where the dynamics is stopped or modified after a given event occurs. 
Concrete examples include leaking systems [l[ and metastable states In some cases we would like to understand 
the probability of a failure event as for example in the spreading of epidemics on networks In other situations, 
hitting a predetermined region in phase space could be a desirable event, as for example if that region is the only 
accessible part from which we can gain information Many of these situations are typically modelled by open 
dynamical systems which we consider here. A closed, discrete-time evolution on a state space X can be opened by 
defining a region E C X in which particles can leak out, early work on such systems can be found in [H, If the 
dynamics is chaotic enough, the measure of points remaining in the system after n iterates will decrease exponentially. 
The escape rate, which is the rate of the exponential decrease, is a well defined object, often admitting an interpretation 
as a spectral quantity and thus invariant for a large class of initial densities [7|, |8| . On the contrary, in [9| it was 
observed that varying the position of the hole has a strong effect on the average lifetime of chaotic transients due to 
the complex periodic orbit structure in chaotic maps. This dependence on the position of the hole has generated a 
renovated interest among both physicists and mathematicians; in [lo| the escape rate was shown to have non-trivial 
dependence on the position and non-monotonic dependence on the size of the holes, general results on the asymptotic 
behavior for small holes have been studied in 11- 3] whilst in 15, T^, the addition of noise is investigated. Open 



problems and reviews of this material can be found in [l|, 0, [l3| . Despite the increasing interest in the properties of 
the escape rate, explicit derivations remain a challenging task. Here we introduce a method which greatly simplifies 
this task for a broad class of systems and that can be used in approximating more general open systems. Via this 
method we derive results for both finite-size and asymptotically small holes. To this end, an explicit formula that 
relates the transfer operators for closed and open dynamical systems together with a return-time operator is derived. 

While escape processes are usually interpreted in the literal sense that points arc removed from the system, we 
consider the following equivalent approach: a point that hits the hole does not leave the system, but rather it is 
coupled to a new point with negative mass. The two coupled particles are then evolved with the same dynamics of 
the closed system and their summed contribution to the total mass vanishes, as should be for an escaped particle. 
While being a mere reformulation of the problem, this approach makes apparent the fact that the effect of a leak is 
fully characterized by its dynamics within the closed system. This idea is built on the work of Lind [isj in topological 
dynamics and is reminiscent of the method of images which can be used to solve certain differential equations with 
Dirichlet boundary conditions in electrostatics and transport problems [l9| . We implement our interpretation of 
escape through the transfer operator of the system. 

Consider a closed dynamics T : X X. Given a density of initial conditions po{x), the Perron- Frobenius operator 
C associated with T evolves it to pn{x) = [C"pa]{x) after n iterates. The transfer operator of the corresponding open 
system is defined by 

CopP ^ - xe)p), (1) 



where xe is the indicator function of the hole E. The open system no longer conserves an invariant measure as 
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the leading eigenvalue A of Cop has modulus smaller than one. The corresponding eigenfunction is the so called 
conditionally invariant density and 7 = — log(A) is the escape rate 0]. 

In the Ulam scheme, the dynamics of T are approximated by a Markov chain in N symbols and the set E is chosen to 
be an element of the fc*'' refinement of the corresponding Markov partition In such situations the appropriate 

operator Eq.([T]) is a transfer matrix written in the basis corresponding to all elements of the k^^ refinement of the 
Markov partition. The size of such a matrix thus grows exponentially with k. By showing that the conditionally 
invariant density is piecewise-constant on a much simpler basis, we construct a transfer matrix for the open system 
which grows linearly with k. Note that this Markov setting and its associated symbolic dynamics has a wide set of 
applications including network dynamics and information theory psj where E is defined as a word (or path) in 
the proper alphabet. 

From now on we consider C to be the operator associated with a given level in the Ulam approximation of a closed 
system, / the corresponding Markov partition with partition parts J;, {I ~ 1,..,N) and its j*'* refinement. Fix 
fc > and choose a hole E E l'^. Equivalently E is defined by a fixed word B of length \B\ = fc + 1. Let po be a linear 
combination of constant functions on the elements of l'^. From Eq.(IT]) we have that 

Pi = jO-opPa = Cpo - CxePq- (2) 

Obviously, Cpo is piecewise constant on and CxePo is constant on the set T{E) and zero elsewhere. Furthermore 

P2 = C^Po - C?XePo - {Cpo - CxEPa) , (3) 

where C'^xePo is piecewise constant on T'^{E) and zero elsewhere. More generally it is easy to see that any function 
that is piecewise constant on any is eventually attracted to the space of functions that are linear combinations of 
constant functions on the elements of | Ji , J/y, T(E), T''~^(E)}. Note that even if these sets are not disjoint, 
their indicator functions are linearly independent [l8|. The space they span is left invariant by the open dynamics 
Eq. (IlJ and thus contains the conditionally invariant measure. We therefore restrict the study of the open dynamics 
to such a space. The operator Cop (with \B\ > 2, but note [sl]) simplifies to 
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(4) 



where A is the N x N transition matrix for the closed system in the original partition / (we assume that the closed 
dynamics is chaotic enough so that A is aperiodic). Denoting hy fj.B and pt the measure of the hole and of its parent 
partition element Jt, with respect to the invariant measure p of the closed system A, we have a = ps/pt- Note that 
the minus signs can be interpreted as couplings to points with negative mass. The 1 in the last row is in the j*'' 
column where Jj is the partition part that the hole gets mapped onto after k iterates, i.e. T^{E) = Jj. The are 
the probabilities that points in E return to E at the i*'' iteration. As we will see, an important role is played by the 
weighted correlation polynomial 



bB{z) = 1 



\B\~2 

1=1 



(5) 



which is a weighted version of the one studied by Guibas and Odlyzko[2^ and Lind |18|] in the context of topological 
dynamics. Note that the coefficients in Eq.© and a are easily computed in all practical situations. 

Not only have we reduced the growth with \B\ of the size of the transfer matrix from exponential to linear, the 
specific form of P allows us to explicitly derive the dynamical zeta function C,op{z) = det(l — zP). The logarithm of 
the smallest zero of this determinant gives the escape rate. Starting from the last row of (1 — zP), multiplying by z 
and adding to the row above and continuing this process up to the {N + 1)*'', we expand the determinant along this 
row and obtain 



(6) 
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FIG. 1: (Colour online) Escape rates for cylinders of length three in the full shift on two symbols 1, with bias p. From p = 3 
vertically we have 000, 010, 001 = 100, 101, 111 and Oil = 110. 



where Cd^(^) ^ dct(l — zA) is the dynamical zcta function of the closed system and Cij{z) = (— dct(l — zA)ij 
denotes the minor on row i and column j. We stress that C~i^{z) and Ct.j{z) are functions of the closed system only 
and are thus independent of the choice of the hole. The escape rate depends on the hole through its returns (via 
<j)B{z)) and through its size (via a and \B\). In particular, finite-size holes with equal measure can have different 
escape rates due to different <f>B{z) while different holes with seemingly different dynamics can have the same escape 
rate if their 4>b{z) and size are equal. In the case of a full-shift we have Cci^{z) = 1 — ^ and Ctj{z) = 5tj{l — z) + ntz 
implying that Eq.([6]) reduces further to 

C-pHz)^{l-z)4>B{z) + fiBZ^''K (7) 

where 4>b{z) differs from (f>B{z) only in that the sum is taken to \B\ — 1 in EqjfS]). 

There is a vast literature on the study of dynamical zeta functions (see [25| and references therein) but general 



expressions that are easy to calculate such as Eqs. ([51 [7]) are often not available [26[. As a particular example of the 
usefulness of having an explicit zeta function, consider a full-shift on two symbols 0, 1 with transition probabilities 
1/p and 1/q such that 1/p + 1/q — 1. From Eq.([7]) we can derive a relation between the dynamical zeta functions of 
the hole defined by the word B ~ O'^' and of its 'child' hole Be = Ol'^'l namely 

Ql,,{z) = QUz){l-'-). (8) 

From Eq.(I5]) we can see that for some values of p the two holes have identical escape rates if the dominant zero 
of Cqibii{z) is equal to that of Cojfl|(-^) while a crossover appears when the leading zero of CqIbiiI^) comes from the 
(1 — z/p) term, as noted in The crossover at p = [i?|/(|i3| — 1) is then obtained by solving Cgj^i (p) = 0- Similar 
results are obtained using the symmetry p— !>(1— p),0— s-l, see Fig. [1] Note that the ordering of the escape rates, for 
finite-size holes, changes with p. At p = 2 (corresponding to the unbiased case) the ordering is given by the length of 
the shortest periodic orbit in the holes (loj . 

An additional strength of Eq.(IB]) is that it allows us to study the asymptotic behavior of the escape rate as fiB shrinks 
to zero around a given point. Consider a periodic point x G A of prime period p (the aperiodic case corresponds 
to p — > oo). Pick a sequence of shrinking holes defined by words Bm chosen to be m concatenated copies of a word 
w — wiW2...Wp. The Wi correspond to the partition sets J^^ visited by the prime periodic orbit of x. Denote by p.^ 
the measure of the interval defined by w. The correlation polynomial of each Bm is then 

ni—l \B \ 1 1 A \ 

(^) - 1 + E = — y^-fjk: — (9) 

n=l ' ^ 

where we use c^p = c^, HBr,^ = c^T^l^w s-nd define A^^ = Cp. For ID maps corresponds to the stability of the 
prime periodic orbit of x. In order to investigate the asymptotic behavior around fiB = e <C 1 we rewrite the smallest 
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zero zq of Eq. (j6l) in terms of a formal expansion 

oo 

0o = l + ^5fee'=. (10) 



k=l 



Using Eq.® and Ea. ((TO| in Eq.® along with the faet that Cd^i^) = (1 ~ z)G{z) for some polynomial G{z) from the 
Perron-Frobenius theorem, to first order in e we have, 



Which gives. 



C(^o) « + -^Cil)t,, = 0. (11) 



In order to make sense of Eq. (|12p we firstly derive an important relationship between the transfer operators of the 
closed and open systems and the operator 

oo 

T{z) - 1 + ^ z"T„ (13) 

n=l 

where r„ is the return time operator of a set E which is defined by Tnp = C^xeP similar to the operators defined in 
[H] . From Eq. ([T]) we have formally 

1 - zLop = 1 - zC + zCxe 

= il-zC)Tiz). (14) 

We stress that this relation is generic and not restricted to the Markov setting investigated so far. On the other hand, 
in the present Markov one, we can show that det(T(z)) = U{z) is the generating function of the probability of return 
to the hole, while det(l — zCop) = Copiz) and det(l — zjC) ~ (^~i^{z). Finally we obtain 



CoV(^) 



U{z) = (15) 



U{z) is related to the generating function of the probability of first-return F{z) via [19|, 

. ^ U(z) - 1 

Using Eq. p^ and Eq. p^ we explicitly show that the decay rate of the first return time distribution of a set _E in a 
closed system, given by the pole of F(z), is equal to the escape rate of the system open on E (given by the zero of 
Copiz)), a previously heuristically derived result Furthermore, from Kac's lemma F'{1) ~ ^/ Pb and thus, using 
Eq. ([T5l) in Eq. (fT6|) we have 



G(l) 



(17) 



which implies 

si-l-A-i. (18) 

Note that for p — > oo (corresponding to x aperiodic) we have A^^ — > and si = 1 (equivalently choose to be 
aperiodic so that t/)B^^(z) = 1 for every m and the result follows). 
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The first order expansion Eq. (|T8|) is consistent with [ll[, wherein it is shown to be valid for more general chaotic 
maps. On the other hand, using Eq.® and Eq. pT)) one can derive all orders recursively due to the general form of 
Eq.®. Of particular note is that for a.e. x Q X, corresponding to x aperiodic, we have, 

. 1 G'(i) c;^ 

From the Entropy (Shannon-McMillan-Brciman) theorem [2^, we have that for /i-a.e. x € X, \B\ — ln(e)//i for 
i?| — > cxD where h is the metric entropy. That is, in the expansion of zq an ln(e) term appears for all subshifts of 



finite type for /i-a.e. x ^ X. The appearance of an e^ln(e) term is derived for the doubling map in jl3j along with a 
heuristic argument for why it should be found in a more general setting, as we have confirmed here for the general 
case of a subshift of finite type. 

Having an expression for the dynamical zeta function such as Eq.® has allowed us to analytically derive some 
properties of the escape rate for a paradigmatic class of dynamical systems. From a mathematical viewpoint, we 
believe that a formal investigation of the accuracy and convergence properties of the proposed framework within the 
Ulam approximation scheme, could permit the extension of some of the results to more general nonlinear systems. 
From a physical viewpoint, the method introduced could be adapted to approximate open systems in many practical 
situations, such as open billiards 1171 . Note that the formalism is equally valid in situations where some piecewise 
constant observable is introduced |l4| or where an infinite countable partition is used as in intermittent maps [30| . 
Finally we stress that the Eas. (l6ll5p allow us to extend Kac-like relations to higher moments of the probability 
distribution of first-return time [3l|: as an example, for a full shift with (pB^z) = we find 

Similar expressions are obtained in the general case. Extending these results will be the subject of future work. 
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